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Abstract. Commuting is an important property in many cases of investiga- 
tion of properties of operators as well as in various applications, especially in 
quantum physics. Using the observation that the generalized weighted differ- 
ential operator of order k and the weighted Hardy-type operator commute we 
derive a number of new and interesting identities involving some functions of 
mathematical physics. 



1 Introduction 

When studying two operators P, Q of quantum theory it is crucial whether the 
relation 

[P,Q]=0, where [P, Q] = (-ih)- x {PQ - QP), 

is fulfilled or not. Such physical quantities, for which there are uniquely assigned 
operators P, Q, are simultaneously measurable if and only if [P, Q] = 0. Equiv- 
alently, we say that the operators P, Q commute if the commutator [P, Q] = 0, 
i.e. PQ = QP. There are many known operators for which the commutation 
relation is fulfilled, e.g. the class of normal operators. On the other side, non- 
commutative operators are source of some interesting stories in mathematics 
and physics, e.g. the famous Heisenberg uncertainty principle saying that it is 
impossible to know the momentum and position of a particle simultaneously 
(see [S]). 

In this paper we are interested in two operators, namely the Hardy-type 
operator H w defined by 

n w f(x) = —^ [ w(t)f(t)dt, x>o, 

where — oo < a < oo, wi(x) = J x w(t) dt and w, / are real measurable, locally 
integrable functions, and the generalized weighted differential operator D™ of 
order k given by 

u>i (x) d \ k . . / w\{x) d\ k 1 / w\(x) 



\ w(x) dx J \ w{x) dx J \ w(x) 

k = 1,2,.... Note that operators H w are of the interest in various functional 
spaces mainly in connection with Hardy's inequality, cf. [4]. The case a = 
and w = const corresponds to the Hardy's averaging operator (or, Hardy's 
arithmetic mean operator), therefore we propose to say H„, the Hardy-type 
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operator. On the other hand, operators appear in the theory of differential 
equations (see e.g. [5] and [7]). If wi(x)/w(x) — 1, the resulting fc-th order 
differential operator will be simply denoted by 



as usual. A generalization of weighted Hardy's averaging operator is provided 
in [3], where the general mean-type inequality involving such operator is inves- 
tigated. 

The first result of this paper is the observation that the operators and 
H w commute. Using this fact in Section 3 we establish few interesting new 
identities involving some special functions of mathematical physics. 

2 A note on commutation relation 

Now we describe an easy observation on commutativity of operators and 
H u ,. For the sake of brevity let us replace H w h(x) by H(x). Then we have 



Differentiating the equality |T]) and then using integration by parts with w\ (a) = 
we obtain the identity 



n k f(x) = —f(x), fe = i,2,... 




(i) 



H'(x) 



w(x) f* Wl (t)h'(t) dt 
Wi(x) 2 



which may be written as follows 



wi(x) 
w(x) 




or in the terms of weighted differential operator as 




Induction on k gives 




k = 



1,2,... . 



(3) 



Summarizing the above we get 



Theorem 2.1 For each k £ N operators D™ and H w commute, i.e. 



(4) 
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Remark 2.2 The observation that operators D™ and H w commute for each 
k e N may be also verified using the well-known fact that the operators 

A/Or) = f{x) + ^f(x) 
w(x) 

and H ffi are the inverses of each other. 

Remark 2.3 Note that for the weighted Hardy's averaging operator (the case 
a = 0) it is enough to consider the operators 



Cf(x) = e~ x [ X eVWdt, 

J — oo 



and L/(ar) = / (lnu>i(x)) for positive function w. Taking into account the 
relations 

LC = H^L, LDfe = DJf L, 

for each k € N, we have the commutativity of operators C and D k , and hence 
H w and D£. 

Further, let us denote Wk{x) = Wk-i(x) dx, where Wo(x) = w(x), and 

r k (x) = f\ ix) , k = l,2,.... 

w k -i{x)w k+1 (x) 

Using wi (a) = and the relation 

(^) />(*)/(*)<*)=»(«) £Mt)f'(t)dt, 

we get 

1 /' ;r 

r 1 (x)—H w h(x) = — r - tui(t)Di/i(t)dt. 

When we repeat this argument fc-times with W£ + i(x) = we(t) dt, for £ = 
1, 2, . . . , k, we obtain 



r fc (*)^(r*_ 1 (*)^(...r 1 ( a! )^H 1B M*) 
1 f x 

— - / Wfe(i)D fe /i(i)di, fc = 1,2,..., 



(5) 



w k + 

which yields the following 

Theorem 2.4 In the sense of the above notation 

RfeH„, = H U)fc Dfe, (6) 
where the operator H k is given by 

R k v(x) = r k (x) ^j- (r k -i(x) ( . . . r x {x) ^- v{x) 
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Remark 2.5 The operator R& will be called the quasi- differential operator of 
order k because it appears in connection with differential equations with quasi- 
derivatives, see e.g. [5]. Clearly, if r k (x) — 1 for each k G N, then R^ reduces 
to Dfc and it is easy to observe that the usual k-th differential operator T> k 
commutes with H w for each k € N. 

Note that the integral identities and ([5]) may be equivalently written as 
differential identities in the form 

Di (w^x^H^x)) = w(x)T>l" (^^(w^x^hix))] , (7) 



and 



Di (w k+ x(x)K k H w h{x)) = w k {x)T> k y— ^F>i(wi(x)li w h(x)J j , (8) 
respectively. Indeed, the term 

1 f x 

H w h(x) = — — / w(t)h(t) dt 
wi{x) J a 

is equivalent to 

— (wi(x)H w h(x) S j = w{x)h{x). 
Now the identity ^ may be rewritten as 

^(wxix^Hvhix)) = w(x)T>%h(x). (9) 



dx V 
If we substitute 

h(x) = w- 1 (x)^(w 1 (x)[ll w h](x^ (10) 
into ©, we obtain the following differential operator identity 



(w 1 (x)T> k B. w h(x)) = w(x)D)» (l-?L( Wl ( x )H w h(x) 
dx \w(x) dx V 

which corresponds to ([7|) . Similarly, the integral identity ((5|) is equivalent to the 
following differential identity 



— w k+ i(x)r k (x) — rk-i(x) — . . . r x (a;) — H w h(x) 
dx \ dx \ dx \ dx 

= w k (x)T) k h{x), 

which may be rewritten by the use of substitution (110[) and the quasi-differential 
operator Rfc to desired form JHJ). 



3 Applications to special functions 

The obtained relations (j4|) and (J6j> involve functions defined as integrals of upper 
limit. Such functions usually appear in the theory of special functions. Special 
functions, also denoted as functions of mathematical physics, have important 
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applications in many areas of mathematics, science and engineering. Using the 
observation on commuting operators we now present few applications of the 
above obtained relations to some special functions to derive some interesting 
identities and representations for them. As far as we know these identities are 
not known in the available literature, e.g. in pQ and [8]. 

In the following few examples we consider a — (the Hardy's averaging 
operator) and polynomial weights. 

Example 3.1 Considering the functions h(t) = e~* and w(t) = t a_1 for a > 
we get wi(x)/w(x) ~ x/a. Then 

1 f x a 

H^O) = / tf^e-* dt = — 7 (a, x), 



fit^dtJo x" 

where 7(0, x) denotes the (lower) incomplete gamma function (see pQ). For the 
generalized k-ih derivative of H w h(x) we get 

1 f d^ k 



r>%H w h(x) = {x— j (ar« 7 (a,aOJ, fc = l,2, 

On the other hand, the direct calculation of H w T)^h(x) yields 

= ^rE(- 1 ) fe " i+1 ^^ fc - i + 1 ) f 

x a , =i J 
1 fc 



t a+k-i e -t dt 



where S^n,™) is the Stirling number of the second kind (see [T]). By Theo- 
rem ET] we obtain the relation 



/ t \ k -1 k 

(x—j (x~ a 7 (a,x)) = — ^(-l) fc - <+1 S 2 (fc, k - i + 1)7(0 + fc - t + 1, x). 

For k — 1 we get the well known relation 7(0 + l,x) — a^y(a,x) — x a e~ x . 
Similarly, for k = 2 we have 7(0 + 2, x) — a(a + 1)7(0, x) — (a + 1 + x)x a e~ x . 

Remark 3.2 In the literature (see pQ) there exists the following representation 
7 (a + n, x) = (-l)«x a +" — (x^ 7 (a, x)). (11) 



Therefore, 

k 

J 



^) " (^ a 7(«, *)) = £>(M-i + l)x fe - < + 1 ^ r (x- 7 (a, 



x)). 
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Example 3.3 Let us consider w(t) = t n , n E No, and h(t) = ipo(t), where tpo(t) 
is the digamma function. Using the same method as in the previous example, 
we get 

T)'kH w h(x) = I ^i FZT (^V*,(,)di) , 



and 



H ro D^/i(a;) 



x n+1 (n 



Since (i^r) V'oW = J2 S2{k,i)t l ipi{t), where ipn(t) = ■^■ipo^t) is the polygamma 

i=l 

function, then 

k 



Thus, 

fe 

a \ 

x- 



Tx) " / * ^ d f ^ fc < ') f r+Vl(<) d " 

Finally, for fc = 1 we get 

J\^Mt)cU = x^^J%^t)*t) 

= x n+1 Mx)~(n + l) [ X t n Mt)dt, 



which gives a formula for the (n + l)-th moment of the trigamma function by 
the use of digamma function. 

Example 3.4 Let w(t) = t m , m € N , h(t) = B n (t), n e N , where B n (t) is 
the n-th Bell polynomial (see [l]). Then 

D™H lu /i(x) = J -±^ (*-§)* (-^ £ ^„(*)dt) , 

and 

H w T>fh(x)^ — — - rr [ t m (t4-) BJt)dt. 

Using the explicit formula 

n 

B n (t)=535 2 (n,i)t i 

i=0 

we get 

1 { x m ( d \ k „ _ . 1 



,„ , / i m t-) g n (t)dt= ? — —— a* 

a; m+1 (m + Jq \ at J (m + \) K ^—^ m + i + l 
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Thus, 



Specially, when m = and fc = 1 we get the relation for integral of the n-th 
Bell polynomial 

Jo + l 

Example 3.5 Consider w(t) = 1 and h(t) = (t 2 - if, k e N. Then we have 

Wk(x) = —, and rjk(a;) = fcW = — — , k = l,2,. 

k\ w k ^ 1 (x)w k+1 (x) k 

A direct calculation yields 

d k (\ 



n k n w h(x) = (k + ^- J (f - \f dt 



and 



H Wk n k h(x) = k ^[ - if dt = 2 -^±^ jf t*p k (t) dt, 

where the Rodrigues representation 

of Legendre polynomials P k (x) has been used (see [T]). Then by Theorem 
we get 

■ fc+l jfc /-i rx 

However, direct calculations yield 



1 -j\t*-ifdt=Y J 

X JO ;_n 



2(fc - i) + 1 



2(fe-i) 



and therefore we have 



dx k \xj y ' ) (2(k - i) + l)(k - 2i)\ 

from which follows the identity 

1 ^ (-l) l (t)(2(fc-*))! 



-2i 



2 fc fc! ^ (2(fc-i) + l)(fc-2i)! 
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Example 3.6 As in the previous example, put w(t) = 1 and h(t) = 
Then 

Rt „„ tW . <* + (_|_ £ e-- d«) _(* + !)£ 

where erf(a;) is the Gauss error function. On the other hand 
U Wk B k h(x) = k -±± ft** (±e-A dt 



1 e 1 



dt k \ a /tt 
2(-l) fc k + 1 



/ t fe i/ fe (i)e-* 2 dt, 
Jo 



V^F x k + x 

where H k is a Hermite polynomial (see pQ). Therefore we get the relation 

Moreover, using the MacLaurin series representation 

[) 0F^n!(2n+l) ' 

v n— v ' 

we get 

Example 3.7 Let = t fe e~*, k S N and to(t) = 1. Then 

H*H.fc(*) = (* + 1)^ jf" t*e"« dt) = (* + 1)^ ( 7(fc + X) 



and 

where the Rodrigues representation 



t k e- l L k (t) dt, 



e 



! dx r ' 



for the Laguerre polynomials L k is used. Then 

x k+i d k / 1 y t + l 



f t k e- t L k (t)dt = 
Jo 



k\ dx k \ x 

According to (fTTj) we have 



dx k \ x J dx k \ x 
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and therefore 

A:! da;* 1 \ x / 4^ \i J il dx k ^ 

Since 



dx k+l \ X I x k+i+l Z-^i ' x3 + 1 

3=0 

where T(n, k) are the permutation coefhcients giving number of permutations 
of n things k at a time, then 



/ t fc e - { L fe (t)dt = y; 

JO ,_ n 



(-1) 4+1 /jfe 



i=0 



fe+i 

e -»^T(A. + i ) j-)a, fc +*--' - (k + i)\ 

3=0 



In the next three examples we derive special function relations involving 
a = +oo, resp. a = — oo, and non-polynomial weights. 

Example 3.8 Let a = +oo and w(t) — e~*. Then Wk(t) — e~* for each k € N, 
and therefore rfc(i) = 1 for each fc S N. Choosing = < a_1 , a > 0, according 
to Remark 12.51 we have 

jfe d k / f°° 

where T(a, cc) = T(a)—j(a, x) — f°° e~*t a_1 dt is the (upper) incomplete gamma 
function (see pQ). Since 

^ j -^r(a,x))=e^£rT(a,x), 

X i=0 X 

then we finally get the relation 

k fji 

^2^-rT(a,x) = (a -l) k T(a- k,x), a > k, 

i=0 

where (x) n = x(x — 1) . . . (x — (n — 1)) is the falling factorial. Observe that the 
case k = 1 in the last relation corresponds to the well-known relation 

r(o + 1, x) = aT(a, x) + x a e~ x . 

(t — m 



Example 3.9 Let a = — oo. Consider functions w(t) — e* and f(t) = exp ^ — 

whcre — 00 < m < 00, and < a < 00. Then by Remark 12.51 we have 



Note that the integral e 2^2 di is the normal distribution function with 
mean m and standard deviation a. Since 

e 2a'i dt = a J - 1 + erf 



2 V V Vlo- 

see [T], then by Leibnitz rule we have 



Using the relation 



d n+1 2 2 

erf(x) = (-l) n ^H n (x)e- x , n e N , 



(see fj]) where H„ is the Hermite polynomial, yields 
d k 



dx k 



(t-m)-» 

e 



\/2cr 



i=0 



On the other hand, 



iffe -e V V =(-l) fc (t-m + ( 7 2 ) fc a- 2fc e 



and therefore 



00 



e * " e y-^-^J dt = (-l) k a- 2k e- x / (t-m+a 2 ) k e J dt. 

dt k ' 



Summarizing the above computations we obtain 



lir 



ix — m\ 2 a- m f * /fc — IV /- ,_ fc+J - (x — m 



Example 3.10 Let a = -00. If w(t) = ie"' 2 / 2 , then wi(t)/w(t) = -1/t for 
t > 0. Putting = sint/t we get 



/ x .2 \ 2 r x 2 

e^V s intdi = eT~ / i 1_ *e"./ fc (t) dt, 
-00 J Jo 



where jk is the spherical Bessel function of the first kind and the relation 

1 d \ k sin i jk (t) 



t dt t t k 
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was used (see pQ). Since 

/ e 2 suit di = ~i\ — erf — =- — err — =- , 
where i is the imaginary unit, then we have 

Similarly, if we choose hit) = t m+1 j m (t), m is integer, and use the Rayleigh's 
relation (see [T]) 

the we get 

(-is) 

4 Concluding remarks 

As far as authors know many of the above obtained identities does not ap- 
pear in the available literature and therefore we suppose they are new. By all 
means the method of their acquirement seems to be new and interesting from 
the viewpoint of applicability of some commutation relations for special opera- 
tors. This method also provides a way how to obtain other possible identities 
involving non-elementary functions (not necessarily special ones) useful in pure 
mathematics, theoretical physics as well as applied sciences. 
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